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ABSTRACT 

It is now generally accepted that long-duration gamma ray bursts (GRBs) are due to the collapse 
of massive rotating stars. The precise collapse process itself, however, is not yet fully understood. 
Strong winds, outbursts, and intense ionizing UV radiation from single stars or strongly interacting 
binaries are expected to destroy the molecular cloud cores that give birth to them and create highly 
complex circumburst environments for the explosion. Such environments might imprint features on 
GRB light curves that uniquely identify the nature of the progenitor and its collapse. We have per- 
formed numerical simulations of realistic environments for a variety of long-duration GRB progenitors 
with ZEUS-MP, and have developed an analytical method for calculating GRB light curves in these 
profiles. Though a full, three-dimensional, relativistic magnetohydrodynamical computational model 
is required to precisely describe the light curve from a GRB in complex environments, our method can 
provide a qualitative understanding of these phenomena. We find that, in the context of the standard 
afterglow model, massive shells around GRBs produce strong signatures in their light curves, and that 
this can distinguish them from those occurring in uniform media or steady winds. These features can 
constrain the mass of the shell and the properties of the wind before and after the ejection. Moreover, 
the interaction of the GRB with the circumburst shell is seen to produce features that are consistent 
with observed X-ray flares that are often attributed to delayed energy injection by the central engine. 
Our algorithm for computing light curves is also applicable to GRBs in a variety of environments such 
as those in high-redshift cosmological halos or protogalaxies, both of which will soon be targets of 
future surveys such as JANUS or Lobster. 

Subject headings: cosmology: theory — galaxies: star clusters — gamma rays: bursts — 
supernovae: general — ISM: clouds — stars:winds, outflows 
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1. INTRODUCTION 

It has long been believed, and pre-explosion progen- 
itors have proven for nearly a dozen cases, that type 
Ib/c and II supernovae are produced in the collapse 
of massive stars. Although we have yet to observe a 
pre-explosion gamma-ray burst (GRB) progenitor, the 
evidence that many long-duration GRBs are also pro- 
duced via massive star collapse has grown since their 
disco v ery in 1973 (for reviews, see IWooslev fc Bloom] 
120061: iFrver et all l2007| ). The engine at the heart of 
these long-duration GRBs is believed to be powered 
either by a rapidly- accreting black hole or a rapidly- 
spinn i ng magnetar (iWoosleyi 119931 IWooslev fe Blooml 
[200l iBarkovfe Komissarovi boTTir but the exact pro- 
genitor of these engines is still unknown. More than 
a dozen scenarios have been proposed in which a mas- 
sive star or tightly- coupled binary system can collapse 
to form such disks (IFrver fe Wooslevl 119981 : IFrver et all 
[19991 [20071: IZhang fc Frverll2001lL 

The simplest progenitor is the classic collapsar, a mas- 
sive star that sheds its hydrogen envelope by strong 
winds and violent mass eruptions (akin to luminous 
blue variable ou tbursts) prior to the collapse of its core 
(| Wooslevl [1993] ) . Additionally, a host of massive star 
models exist, a set of which invoke binary merger events. 
Like the violent mass ejection in massive stars, these 
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merger events eject shells of material into the immediate 
surroundings before the GRB outburst. Both single star 
and binary mass ejections can occur just a few thousand 
years prior to stellar collapse. In the case of the helium 
merger model, we expect the merger-driven mass ejec- 
tion to occur less than a few years before the launch of 
the GRB jet. The jet must plow through this shell as it 
is producing the gamma-ray emission we observe. 

The structure of the circumburst media (CSM) should 
imprint signatures on their afterg l ows that identi fy the 
mode of collapse ([Meszarosl l2002t IWooslev! 120 111) , pro- 
viding a mechanism that can be used to better un- 
derstand the progenitors of these cosmic explosions. 
Analytical models of relativistic jets in both uniform 
circumburst densities and free-streaming wind profiles 
have yielded light curves that are in reasonable agree- 
ment with observations (iBerger et all 120001 : lYost et all 
12001 iCurran et all [201lt iPrice et al.M2002l ). Recently, 
collisions of the jet with more complicated struc- 
tures h ave been examined, such as wind-termination 
shocks (iRamirez-Ruiz eFaIll200iT 12001 iDai fe LvJ 120021: 
INakar fc Granotll2007lh shocks due to collisions between 
stella r winds f rom nearby stars. dMimica fc Gianniosl 
[20T1 . c lumps (E amirez-Ruiz et al. 1I2005IL and magnetic 
shock s (I Yost et all l2003lh IRamirez-Ruiz et al.1 (|200lL 
I2005| ) and IDai fc Lul (|2002l ) found that sharp features 
in circumburst densities due to clumps and shocks cre- 
ate discernible features in GRB light curves, although 
INakar fc Granotl (|2007ft . who consider the dynamics of 
the reverse shock and assume that the jet remains rel- 
ativistic after encountering a jump in density, do not. 
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iMimica fc Gianniosl (|2011[) found gradual shallowing in The hydrodynamics equations are 
the light curve when the jet passes through a shock 
formed when stellar winds from the progenitor and a 
neighbor collide. 

Flares are often observed in GRB lightcurves within 
the first few thousa nd seconds, nota bly in GRB 08102 9 
(INardini et al.ll20rl. GRB 0711 12C (iHuang et al.ll20"H , 
GRB 050502B (iFalconeeFall 12006]), G RB~060607A 
tiZiaeepour eTall 120081) . and GRB 050421 (jGodet et al.l 
120061 ). It is generally accepted that these flares are a 
product of late-time injections of energy into the system 
by the GRB's central engine, but this is not the only 
possible scenario. When the jet from a GRB encounters 
a shell produced by a stellar progenitor, it experiences 
an abrupt increase in the medium density. GRB light 
curves are much more sensiti ve to a decrease in medium 
density than to an increase (jYost et al.l 12001 . meaning 
that an abrupt increase in density of the order of ~ a few 
will leave almost no discernible imprint on the observed 
light curve. In contrast, our hydrodynamical models can 
produce sudden enhancements in the density of five or- 
ders of magnitude or more at the trailing edge of the 
shell (Fig.Qj, which is sufficient to produce bright flares 
despite the weak dependence of the light curve on an 
increase in density. 

Winds, outbursts and ionizing UV radiation from GRB 
progenitors all disperse the molecular cloud cores that 
created them and create far more complex ambient mor- 
phologies for the jet tha n those considered in afterglow 
studi es to date (but note lFrver et al.l[2006l : lw"halen et al.1 
[2008]) . In this paper, we model such environments for a 
variety of GRB scenarios with ZEUS-MP. Our ID models 
span a wide variety of winds and outbursts with complex 
gas chemistry and cooling that capture the true structure 
of the circumburst medium. We have also developed 



-V • (pv) (1) 
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where p, e, and the Vi are the gas density, internal energy 
density, and velocity of each zone and p = (7 — 1) e and Q 
are the gas pressure and the von Neumann-Richtmeyer 
artificial viscosity tensor. We evolve mass fractions for 
H, H+, He, He+, Hc 2 +, H , H+, H 2 , and e" with nine 
additional continuity equations (the species are assumed 
to have the same velocities) and t he nonequilibrium rate 
equations of lAnninos et al.l (|1997t) 



^ = -V • (pv) + £ Pjk(T) PjPk + £ K jPj , (4) 



where f3j k is the rate coefficient for the reaction be- 
tween species j and k that creates (+) or destroys (-) 
species i, and the the radiative reaction rates. 

Microphysical cooling and heating are calculated with 
operator-split isochoric updates to the gas energy den- 
sity that are performed every time the reaction network 
is solved: 



A. 



(5) 



a semi-analytical app roach b ased on previ ous work by 
Panaitescu fc Kumar! (|2000l) . iHuang et all (|1999D . and 
Pe'erl (|2012l ) for computing GRB light curves in any gen- 



eral density profile, not just the uniform media, winds 
and simple density jumps in previous studies. We apply 
this new method to compute light curves for rclativistic 
jets propagating through circumburst shells, and exam- 
ine the imprint of these shells on the light curves in order 
to determine if they constrain the mode of collapse. 

In § 2 we discuss the ZEUS-MP code, how it is used to 
simulate the environments of long-duration GRBs, and 
our grid of wind models. We review the results of our 
shell ejection calculations in §3. In §4 we describe how 
our new analytical models of GRB jets are applied to 
hydrodynamical profiles from these simulations to com- 
pute light curves for a variety of energies. We calculate 
GRB light curves in our grid of shell profiles, which cor- 
respond to a variety of collapse scenarios, and determine 
if specific light curve features constrain the nature of the 
progenitor. In § 5 we conclude. 

2. NUMERICAL METHOD 

2.1. ZEUS-MP 

ZEUS-MP is a massively-parallel astrophysical hydro- 
dynamics code that solves nonequilibrium H and He gas 
chemistry and photon-conserving ionizing UV radiation 
transport together with Eulerian fluid dynami c s in a s elf- 
consistent manner (Whale n fe Norman! 120061 l2008al lbh . 



Here, T is the photoionization heating rate for all 
species over all photon energies and A is the sum of all 
cooling rates. We include collisional excitation and ion- 
ization cooling by H and He, recombinational cooling, H2 
cooling, and brcmsstrahlung cooling, with our nonequi- 
librium reaction network providing the species mass frac- 
tions needed to accurately calculate these collisional cool- 
ing processes. We also calculate fine structure cooling 
due to C, O, N, Si and Fe using the lDalgarno fc McCravl 
(|1972h cooling curves, generalized to arbitrary elemental 
abundances. We exclude cooling by dust. 

Fluid flow, gas heating and cooling, and H and He 
chemistry can occur on highly disparate timescales whose 
relative magnitudes can widely vary throughout the 
course of a calculation. The many chemical reaction 
timescales can be consolidated into a single chemistry 
time step defined as 



t 



chem 



= 0.1 



(6) 



which is formulated to ensure that the fastest reaction 
operating at any place or time on the grid determines 
the maximum time by which the reaction network may 
be accurately advanced. The timescale on which the gas 
heats or cools is given by 



0.1 



(7) 
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To evolve each physical process on its respective 
timescalc without restricting the entire algorithm to the 
shortest one, we subcycle the reaction network and en- 
ergy equation over the time step on which we update the 
hydrodynamics equations. We first compute the mini- 
mum t c hem and th/ c for the entire grid and then perform 
consecutive updates of species mass fractions and gas en- 
ergy densities over the smaller of these two times until 
the lesser of th/ c an d tcFL is reached, where tcFL is the 
Courant time. At this point a full update of the hydro- 
dynamics equations is performed and the cycle repeats. 

The prcfactor of 0.1 in equations [5] and [7] guarantees 
that mass fractions never change by more than 10% when 
the reaction network is solved and that gas energies do 
not change by more than 10% over a time step. This pre- 
vents catastrophic runaway cooling in gas at high densi- 
ties and metallicities by either H2 or fine-structure cool- 
ing, like those in dense shells swept up by strong winds 
or mass ejections. 

2.2. Problem Setup 

To model mass ejections in ZEUS-MP, we treat stellar 
winds and outbursts as time-dependent inflows at the in- 
ner boundary of a ID spherical grid with 32,000 zones. 
The gas is assigned H and He mass fractions of 0.76 and 
0.24, respectively, and a metallicity Z = 0.1 Z@. The 
mesh extends from 10~ 4 pc to 0.3 pc with outflow con- 
ditions on the outer boundary. The inflow is imposed at 
the inner boundary in the form of a time- varying density 
and velocity: 



P 



4:Trr lb 2 v w ' 



(8) 



where is the radius of the inner boundary and v w 
is the wind velocity. Outbursts are modeled by increas- 
ing m and lowering v w . Because stellar winds clear out 
gas from the vicinity of the star prior to any outbursts, 
we initialize the grid with a free-streaming density and 
velocity profile 



p{r) 



Airr 2 v v 



(9) 



where the wind velocity is assumed to be constant. 
The temperatures of the initial density profile is set to 
100 K. We launch the outburst at the beginning of the 
simulation. The grid is domain decomposed into 8 tiles, 
with 4000 mesh zones per tile and one tile per processor. 

We neglect the effect of ionizing radiation from the 
star on the structure of the dense shell. This treatment 
is approximate, given that the progenitor illuminates the 
flow over its entire lifetime and that its luminosity evolves 
over this period. However, the heat deposited in the wind 
by photoionizations is small in comparison to its bulk 
kinetic energy and is unlikely to alter the properties of 
the flow in the proximity of the GRB. 

2.3. Grid of Shell Models 

For collapsar and He mergers we consider 3 mass loss 
10~ 6 ,10- 5 , and 10~ 4 M Q /yr and outbursts 



rates m„ 



10~ 2 M Q /yr lasting for 10 yr and 100 yr that 



correspond to total shell masses of 0.1 and 1.0 M Q , re- 
spectively. We take the velocities of the fast wind and 
slow shell to be 2000 km/s and 200 km/s, respectively, 
and in each model use the given rh w to initialize the den- 
sity across the entire grid, assuming that the star sheds 
mass at the same rate before and after the outburst. In 
general, larger shell masses are expected for He mergers 
but both kinds of progenitors can exhibit light, moder- 
ate and heavy winds, so there is some degeneracy across 
our grid of models. We also consider He-He mergers by 
simulating the loss of the common envelope with a single 
massive outburst rn^ = 10 M Q /yr that lasts for one year 
and has a velocity of 200 km/s. 

2.4. Wind Bubble Test 

Before calculating circumburst environments for a 
GRB progenitor, we first model the bubble blown by its 
fast wind during the life of the star. This bubble is the 
primary circumstellar structure formed by the star and 
grows to radii of 20 - 30 pc before outbursts alter its in- 
terior on scales of a few tenths of a parsec late in the life 
of the star. Since the gamma-ray spectra and afterglow 
of the GRB are primarily governed by the interaction of 
the jet with its surroundings out to ~ 0.1 pc and 10 pc, 
respectively, this primary shell does not affect the obser- 
vational signature of the GRB at early times. However, 
we simulate this bubble first to verify that for reasonable 
choices of ambient density it is indeed at least 10 pc from 
the star when it dies and does not immediately impact 
the afterglow. Also, the effect of local density and ra- 
diative cooling on the structure and kinematics of shells 
in general is more easily seen with this first bubble than 
the more complicated structures created by interactions 
between slow outbursts and fast winds just before the 
death of the star. 

The wind bub ble has the cl a ssic t wo-s hock structure 
first d escribed bv lCastor et al.l ()1975h and lWeaver et al.l 
(fl97l . As we show in Fig. [1] one shock forms at the 
interface between the emergent wind and the surround- 
ing gas as it is swept up at speeds greatly exceeding the 
sound speed of the gas. As gas accumulates on the bub- 
ble, the shock detaches from the wind and moves ahead 
of it, forming an intervening shell of dense postshock gas. 
At the same time, the expansion of the bubble evacuates 
a cavity into which the wind freely streams. Since the 
shell moves more slowly than the wind, a termination 
shock also forms where the wind piles up against the in- 
ner surface of the shell. If gas in the shell can radiatively 
cool, it flattens into a cold, dense structure that is prone 
to fragmentation into clumps. 

We first performed four tests with steady winds of rh w 
= 10~ 5 M Q /yr and v w = 1000 km/s in uniform densities 
n = 10, 100, 1000 and 1.8 x 10 4 cm" 3 to investigate 
how local densities govern the radius of the bubble at in- 
termediate times. For simplicity, these calculations were 
done with no chemistry or radiative cooling. As we show 
in the left panel of Fig. [I] ambient density governs only 
how far the bubble is driven from the star, not the profile 
of the free-streaming region in the immediate vicinity of 
the star, which is determined only by m and v w (and 
fluctuations thereof). 

We find that radiative cooling has a dramatic effect on 
the structure of the primary shell but no influence on the 
flow up to the shell, as we show in the right panel of Fig. [I] 
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for the same wind and n = 100 cm" 3 . The three plots 
show the structure of the bubble with no cooling, H2 
cooling, and fine-structure cooling due to metals at Z = 
0.1 Zq. The cooling flattens the shell into a cold dense 
structure and radiates away some of its thermal energy, 
slightly retarding its advance. The free-streaming zone is 
again unaffected because the wind velocity ensures that 
the termination shock is at least 0.6 pc from the star by 
the end of the simulation. Thus, chemistry and cooling 
will clearly cause any subsidiary shells ejected by the pro- 
genitor at late times to be much thinner and denser, with 
potentially important consequences for the propagation 
of the jet. 

3. CIRCUMBURST DENSITY PROFILES OF COLLAPSARS 
AND HE MERGERS 

We now consider the more complicated structures in 
the immediate vicinity of the star at the time of the GRB. 
We examine a fiducial case from our grid of shell models, 
the 100 yr outburst in a stellar wind with m w = 10~ 5 
MQ/yr, whose density and temperature profiles we show 
at 120, 600 and 1000 yr in Fig. [2] As the shell emerges 
from the star, it promptly cools to extremely low tem- 
peratures, as we show in panel B of Fig. [2] at 120 yr. 
This happens because fine-structure cooling timcscales 
are less than a year in the dense 0.1 Zq shell. At the 
same time, the fast wind detaches from and races ahead 
of the shell, as we show in panel A of Fig. [2l This creates 
a zone of rarefaction into which the shell freely streams, 
as shown by its roughly r~ 2 density profile. 

As the wind pulls away from the shell, the abrupt 
adiabatic expansion causes temperatures at its inner 
edge to drop sharply. The expansion of the rarefied re- 
gion likewise causes its temperatures to fall to ~ 1 K. 
In reality, gas temperatures cannot fall below the cos- 
mic microwave background temperature at this epoch, 
Tcmb = 2.73(1 + z) K, but this is not included in our 
simulations. We also ignore any heating of the gas due 
to other external processes such as cosmic rays or pho- 
tons from the progenitor itself. Although these processes 
could conceivably alter the density profile in the regions 
of free-streaming wind if those regions were to become 
ionized, the temperature of the shell is set by the conver- 
sion of the free-streaming wind's kinetic energy into ther- 
mal energy at the shell's inner edge (as well as through 
radiative cooling throughout its interior). Any modifi- 
cations of the shell structure due to heating by other 
processes is negligible. 

The density dip at the beginning of the rarefaction 
region and the density jump at the rear edge of the de- 
tached wind arc complementary and due to mass conser- 
vation. When the fast wind breaks away from the shell 
a thin shell of dense gas from its outer layers breaks off 
with it, leaving a thin layer whose density is even lower 
than that of the rarefied region and which remains a per- 
sistent feature of the flow out to ~ 500 yr. A transient 
structure with multiple strong density jumps results: the 
emerging massive shell and the inner edge of the rapidly 
receding wind separated by an intervening low density 
region. The shell soon begins to plow up the low-density 
gas, as can be seen in the density bump just beyond its 
leading edge. 

After the shell has been fully ejected (and the fast wind 
has exited the grid), it soon evolves into a wind bubble 



as the fast wind piling up at its inner surface forms a ter- 
mination shock that detaches and recedes from it in the 
frame of the shell. The termination shock heats the gas 
piling up at the inner surface of the shell to ~ 10 8 K while 
the shell itself remains cold due to fine-structure cool- 
ing. Densities in the postshock gas between the termina- 
tion shock and the shell become roughly uniform because 
sound speeds in the shell are ~ 1100 km/s and it is nearly 
isothermal, so pressure gradients in this region arising 
from any initial density gradients across it are erased by 
acoustic waves on timescales that are short in compari- 
son to the expansion times of the shell. Temperatures in 
the shell are lower than those of canonical wind bubbles 
because the ambient density is too low to heat the shell 
and activate H and He line cooling, in contrast to Fig. [TJ 
Strong brcmsstrahlung x-ray flux from shocked gas at 
the inner surface of the shell likely ionizes it to some 
degree, but probably does not otherwise alter its proper- 
ties because, at these photon energies, most of the kinetic 
energy of the photoelectro ns goes into secondary ioniza- 
tions rather than h eat (e . g. iShull fc van Steenberdll985l ; 
iRicotti et al1l200ll I2002T |2005j). By 760 yr the shell has 
swept up enough low density gas to form a shock, which 
soon separates from the shell and advances beyond it as 
shown in the density and temperature plots at 1000 yr. 
This feature is temporary: the shell eventually subsumes 
the secondary shock as it expands. A free-streaming re- 
gion forms behind the termination shock and extends to 
~ 0.1 pc by 1000 yr. Part of the reason the shell drives a 
shock at intermediate times is that it accelerates as it is 
pushed by the fast wind and as it expands into low den- 
sities that continue to fall over time, as we show in all 
three density plots. A velocity gradient develops across 
the shell as its inner and outer surfaces accelerate to 220 
km/s and 280 km/s respectively over 1000 yr. 

Radiative cooling flattens the ejected shell as it ex- 
pands into the surrounding medium. Its width at first 
is 0.02 pc at 100 yr but later decreases essentially to the 
resolution limit of the grid by 1000 yr. This is not un- 
expected: when a dense shell driven by winds or swept 
up by a supernova can radiativcly cool, its width usually 
drops to a few mesh zones, with the number of zones be- 
ing determined by the numerical viscosity of the hydro- 
dynamic algorithm. All seven models exhibit the same 
general evolution, with the only variations being in the 
thickness of the shell over time, its peak densities, and 
its location on the grid at the time of the GRB. De- 
pending on when the explosion occurs, the jet will en- 
counter three density jumps or four, with jump ratios 
of up to 10 10 . The thickness and average density of the 
shell varies from 0.02 to 0.0001 pc and 10 5 - 10 9 cm" 3 
over 1000 yr. 

All our models of massive shell ejection follow the evo- 
lutionary sequence we have just described, with the only 
differences being in the magnitudes of density jumps, the 
mass and thickness of the shell at a given time, and the 
positions of the regions on the grid. 

4. AFTERGLOW LIGHT CURVES 

We now discuss the general method we have devised for 
computing light curves for relativistic jets in the compli- 
cated wind structures we have modeled in ZEUS-MP. 

4.1. Blast Wave Hydrodynamics 
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Fig. 1. — Left panel: wind-blown bubbles at 6.25 X 10 4 yr for m w = 10 5 MQ/yr and v w = 10 3 km/s for four fiducial ambient densities, 
10, 100, 100 and 1.8 X 10 4 cm" 3 . These profiles include no radiative cooling. As can be seen in the plots, densities within ~ 1 pc of the 
star depend only on rh and v w for n < 100 cm~ 3 . Right panel: structure of the shell at 1.25 X 10 4 yr in an ambient density of 100 cm -3 
with no cooling, H2 cooling and fine structure cooling by C, O, N, Si and Fe at Z = 0.1 Z@ (from the Dalgarno-McCray, or DM, cooling 
curves). Radiative cooling flattens the shell plowed up by the wind into a cold dense structure, with no effect on the free-streaming region 
in the vicinity of the star. Note also that efficient cooling in the shell also radiates away some of its thermal energy and slows its advance. 
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Fig. 2. — Density profiles (a) and temperature profiles (b) for 



a 100 yr outburst with rhj, 



10" 



M©/y r in a stellar wind with 



10" 



M /yr. Black: 120 yr; red: 600 yr; blue: 1000 yr. 



In the canonical fireball model, gamma ray bursts are 
modelled as initially highly-relativistic, adiabatic jets 
that propagate outward into an ambient medium. Wc 
will assume that the jet expands adiabatically, i.e., only 
a very small fraction of the total burst energy is available 
to the electrons to be radiated away, and that all of the 
circumburst material in the jet's path is swept up. In 
the following discussion, primed quantities refer to the 
reference frame that is comoving with the jet, unprimcd 
quantities with no subscript refer to the reference frame 
in which the ISM is at rest, and quantities with the sub- 
script refer to the reference frame of an Earthbound 
observer. 

Energy conservation requirements yield a formula for 
the evolution of the jet as it propagates through the ex- 
ternal medium. In the case where the jet expands adia- 



batically, 

dT 

dm 



7(r 2 -i)-(7-i)r/3 2 

Af cj +m[2 7 r-( 7 -l)(l + r- 2 )]' 



(10) 



where T is the Lorentz factor of the jet, M e j is the initial 
mass of the jet ejecta, m is the total mass that has been 

1 /2 

swept up by the jet, /3 = (l — T~ 2 ) in the normalized 
bulk velocity, and 7 is the adiabatic index (jPe'erll20"T2h . 

The high resolution of our simulations (10 -4 pc) allows 
us to take the density to be constant across each mesh 
zone. The total mass swept up by the jet by the time it 
reaches grid point n is then approximately 



M{r) = -7T 



Pir\ 



j=2 



(11) 



where and /?, are the radius and density of the ith 
grid point, respectively. The time t \, s at which a photon 
emitted at the shock boundary reaches an observer along 
the line of sight can b e calcu lated by integrating equation 
12 from lHuang et all (|1999f) : 



1 

t = - 



dr 



Vr 2 



1 



(12) 



where j3 = v/c is the jet velocity and c is the speed of 
light. 

Many GRBs are thought to produce collimated jets 
rather than isotropic outflows. The center of a relativistic 
jet of half-opening angle 9j is not in causal contact with 



its edge until < D b s 



jet - 



which is defined as the time 



where T ~ 1/%- The jet evolves in an identical manner 
as in the isotropic case until i b s = £jct, at which point 
the jet experiences rapid lateral expansion. The increase 
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in external mass being swept up by the jet after i b s = 
tjet causes the jet to decelerate at an increasing rate. 
Additionally, the afterglow's decreasing luminosity is no 
longer being partially offset by an increase in the size of 
the emitting region as seen by the observer, leading to 
a break in the light curve. The overall evolution of the 
jet's angular size is 



alt 



(13) 



where c' s is th e comoving sound speed and r is the 
radius of the jet (|Huang et al.ll2000D . 

4.2. The Injection Break 

If we assume that a constant fraction eg of the to- 
tal fireball energy is stored in magnetic fields, then the 
equipartiti on magnetic field strength at the shock bound- 
ary is (i. e. iPanaitescu fc Kumarll2000D : 

^ = Ae B m p c 2 n(r)(r - 1) (t + , (14) 

where m p is the proton mass and n(r) is the number 
density of the medium at radius r. 

The electrons that are injected into the shock are as- 
sumed to have a velocity distribution AT (7) oc 7~ p with 
a minimum Lorcntz factor 7 m . Electrons with a Lorentz 
factor 7 e e mit synchrotron radiation at a characteristic 
frequency (|Rvbicki fc Lightmanlll979l ): 



D 



2irm P 



(15) 



7n 



The injection break, is m , corresponds to the characteris- 
tic frequency at which the electrons having the minimum 
Lorentz factor ra diate. The minimum Lorentz factor is 
(ISari et al.lll998h 

f^)^ee(r-l), (16) 
\p - I J m e 

where m e is the electron masses, respectively. Note that 
equation [16] is only valid in the rclativistic limit. Sub- 
stituting equation [TB] into equation ll5 | will therefor e only 
yield v m as long as V > 2.0. iFrail et al.l (|2000D show 
that, for r < 2.0, the injection break frequency follows 
the simple relation v m oc i~ 3 . 

4.3. The Cooling Break 

Rclativistic electrons in the shock cool radiatively 
through inverse Compton (IC) scattering and syn- 
chrotron emission on a co-moving frame timescale 



4d(7) 



6tt 



in,,c 



Y + 1 a e jB' 2 ' 



(17) 



where Y is the Compton par ameter and o~ e is the Thomp - 
son scattering cross section (jPanaitescu fe Kumarll2000h . 
An electron with Lorentz factor 7 C cools radiatively on a 
timescale equal to the current age of the remnant. Solv- 
ing Eqn. [T7]for 7 C , we find that the Lorentz factor for 
electrons that cool on a timescale equal to the observer- 
frame age of the remnant is 



67TTO e C 2 



B' 2 a e (Y + l)t' 



(18) 



4.3.1. Fast-Cooling Electrons 



Electrons in the GRB jet can cool by adiabatic expan- 
sion of the gas or by emission of radiation. When the 
cooling timescale for electrons with Lorentz factor j m is 
less than the age of the jet (y c < v m , where v c is the fre- 
quency of the cooling break) the electrons in the jet lose 
a significant portion of their energy through emission of 
radiation and are said to be radiative, or fast-cooling. 
Conversely, if the cooling timescale is greater than the 
age of the jet (v c > v m ) the electrons do not lose signifi- 
cant energy to radiation and are said to be adiabatic, or 
slow-cooling. 

To calculate the Compton parameter, Y, we only ac- 
count for one upscattcring of the synchrotron photons. If 
the injected electrons are fast-cooling and the frequency 
of the absorption break v a < min(t/ m , v r ), then Y can be 



approximated bv IPanaitescu fc: Meszarosl (j2000f ): 



Y r = 7m7cTc, 



(19) 



where a constant of order unity has been ignored and r e 
is the optical depth to electron scattering, given by 



Q-eMjr) 
47tm p r' 2 



(20) 



The medium becomes optically thick to synchrotron self- 
absorption at the absorption break frequency v a . When 
both the injection break and the cooling break lie in the 
optically thick regime, Y becomes 

Y r = Y, = ri(ctH^- 1) ) 1/(P+5 \ (2D 

where C2 = hq e T' e /o- e B' (jPanaitescu fc Meszarosl[2000t ). 

4.3.2. Slow-Cooling Electrons 
If the electrons arc slow-cooling, then Y becomes 



Y, 



Tell 



(22) 

as lo ng as v a < mm(v m , v c ) (jPanaitescu fc Meszarosl 
2000) and 1 < p < 3. Once again we have ignored a 
constant of order unity. If both the injection and cooling 
breaks lie in the part of the spectrum that is optically 
thick to synchrotron self-absorption, then Y is identical 
to the corresponding fast-cooling case and 



Y, = Y, 



4.4. The Absorption Break 



(23) 



At lower frequencies, the medium through which the 
jet propagates becomes optically thick to synchrotron 
self-absorption. The result is a transition to a F v oc v 2 
drop-off in the flux at some absorption break frequency 
v a where the optical depth to self-absorption is r a b = 1. 
The frequency of the absorption break depends on the 
electron cooling regime (fast or slow) and on the order 
and values of both the i njection and cooling bre a ks. In 
the fast-cooling regime, IPanaitescu fc Meszarosl (|2000D 
find that 
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(-3/10 -1/2 



l/(p+5) 



7a < 7c < 7n 
7c < 7a < 7n 



(C 2 7c7rn _1 )" > 7c < 7m < 7a, 

(24) 

whereas in the case where the electrons are slow-cooling 



f £.3/10^ -1/2 



slow- 1 



7n 



l(C 27 ^ 1 7c) 



1/(p-4) 



l/(p+5) 



7a < 7m < 7c 
7m < 7a < 7c 
7m < 7c < 7a- 

(25) 



4.5. Light Curves 

In order to produce light curves, we must first find the 
time dependence of T(r), n(r), and M(r). Equation HD1 
can be solved numerically for T(r), and equation 1121 can 
then be used to relate the observer time f b s to the jet po- 
sition r, allowing us to rewrite the equations defining the 
three break frequencies in terms of i bs, r(i bs), n(t bs), 
and -M(i b s ). Given the three break frequencies and the 
peak flux density, analytical light curves can then be cal- 
culated that are valid from the radio to the 7-ray regions 
of the spectrum. If v a < imn(v m ,v c ), then the peak flux 
density F™% x occurs at the injection break if v B1 < v c 
and at the cooling break if v m > v c : 



F" 



V34> v 



TB'M{r) 



AttD 2 



m e c* 



(26) 



where 4> p is a factor calculated by iWijers fc Galamal 
(|1999| ) that depends on the value of p, and D = (1 + 
z)~ 1 / 2 Di, where Dj. is the lumin osity distance to the 
source (|Panaitescu fc Kumarll2000| ) . The flux at any fre- 
quency v (ignoring relativistic beaming and the spheri- 
cal nature of th e emi tting region) has been deriy ed by 
iSari et al.l (119981 ) and iPanaitescu fc Kumarl (|2000f ). 

4.5.1. Fast-Cooling Electrons 

When the electrons are in the fast-cooling regime, the 
peak flux density occurs at the cooling break as long as 
v a < v c : 



In the canonical afterglow models that assume a uni- 
form density environment, the cooling break and the in- 
jection break move to lower frequencies with time. Even- 
tually, both the cooling break and the injection break 
can lie below the absorption break, but far too late in 
the evolution of the burst to be relevant to anything but 
the radio afterglow, and long after the time at which the 
electrons in the jet have transitioned to the slow-cooling 
regime. In our more realistic density profile models, the 
extremely high density encountered by the jet as it passes 
through the thick shell causes it to abruptly transition 
from highly relativistic to Newtonian expansion. The de- 
crease in r leads to a sharp drop in the injection break 
frequency, while the increased medium density leads to 
a larger magnetic field strength, which in turn causes a 
drop in the cooling break frequency. The result is that 
the absorption break frequency can be several orders of 
magnitude higher than the cooling and injection break 
frequencies as the jet traverses the thick shell. Multi- 
ple transitions between fast and slow electron cooling 
can also occur. In the vicinity of the thick shell, when 
v a > v m and the electrons are in the fast-cooling regime: 

( (v/v c ) 2 (v c /v a ) 5 ' 2 , v < v c 
F v ^=F^Uvlv a fl\ v c <v<v a (29) 

l(*v^r p/2 , V m < V- 

4.5.2. Slow-Cooling Electrons 

Our models yield the same flux as the canonical wind 
models until the jet encounters the shocked wind that 
has piled up behind the thick ejecta shell. If it does 
not encounter the shocked wind in the first few hours, 
the electrons in the shock transition to the slow-cooling 
regime, with v a -C v m and 

' {v/v a ) 2 {v a /v m Y/ 3 , v<v a 

(tVfm) 1/3 , v a < V < U m 

(tVfm)~ (P " 1)/2 , V m < V < V c 

\{vlv c )-vl\v c lv m )-^l\ V C <V. 

(30) 

In the case where v m < v a < u c , 



F" 



p iTimax 



'{v/v a f(v a /v c ) 1 ^, v<v a 

(i//^) 1 / 3 , V a < V < v c 

(v/v c y 1/2 , V c < V < v m 

Xv/v m )- p/ \Vm/Vc)- 1/2 , V m <V. 

(27) 

If the medium is optically thick to synchrotron self- 
absorption at the cooling break frequency, then the max- 
imum flux moves to the absorption break frequency. 
Between the absorption break and the cooling break, 
F v cx v°l 2 but it becomes cx v 2 below the cooling break: 



F, 



F D 



{v/v c ) 2 { Vc /v a f/ 2 , 
{vlv a fl\ 

{V/Va)- 1 ' 2 , 

(^m)- p/2 K>a)- 1/2 , 



v < v c 
V c < V < v a 
V a < V < v m 
v m < V. 

(28) 



F, 



F*! 



'{v/v m ) 2 {is m /is a ) 5 / 2 , v<v a 
{v/v a f/ 2 , V a < V < v m 

(iV^)-(p- 1 )/ 2 , v m <v<v c 

{(v/v c )-P/ 2 (is c /v a )-(P-V/ 2 , v c <v. 

(31) 

As noted earlier, as the jet passes through the thick 
shell, it can experience multiple transitions between fast 
and slow electron cooling. When the electrons are in the 
slow-cooling regime and v a > v c , 

' {v/v m ) 2 {v m /v a f/ 2 , V <v m 
K {v/v a f/ 2 , v m <v<y a 
>K)- p/2 , v a <v. 



(32) 
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Fig. 3. — Light curves and break frequencies for a GRB in a dense 
shell. Panel (a): densities encountered by the jet over time; panel 
(b): synchrotron light curves; panel (c): break frequencies. Region 
V, a second region of unshocked wind, is not shown here because 
the fast detached wind has exited the grid. 

The spherical nature of GRB afterglow emis- 
sion can substantiall y alter the observed light curve 
(|Fenimore et al.|[l996() . The burst ejecta is initially ultra- 
relativistic, with 100 < T < 1000, meaning that radiation 
that is emitted along the observer's line of sight will be 
beamed toward the observer more than radiation that is 
emitted away from the line of sight. There will also be 
a delay in the arrival time of photons that are emitted 
from regions of the jet that lie away from the observer's 
line of sight, since these regions ar e further away from 
the observer (|Fenimore et al.lll996t) . The overall effect 
will be to increase the total observed flux at early times, 
and to make the light curve broader and more smooth. 

4.7. Imprint of Dense Shells on GRB light curves 

The shell and wind enclosing the GRB can be parti- 
tioned into four distinct regions, which we show in Fig.|3K 
for a progenitor with a 0.1 M Q shell ejected 100 yr before 
the burst and mass loss rates of 10 -6 M /yr before and 
after the ejection of the shell. A fifth region, not shown in 
Fig. [31 also exists, but is located at such large radii at the 
time of the GRB that the jet does not encounter it until 
very late times. Even so, this region is discussed below 
for completeness. We now examine the imprint of each 
region on the GRB light curves and break frequencies, 
which we show in Figures \5jp and|3b. 

4.7.1. Region I - Unshocked Wind 

Region I is the unshocked wind blown by the progenitor 
after the ejection of the shell and prior to the burst. Not 
surprisingly, the light curve in this region is what would 
be expected for a p(r) tx r~ 2 wind profile (Fig. [3J3) . The 
jet electrons are initially in the fast cooling regime, and 



the absorption break frequency is generally higher than 
the cooling break frequency for the first few seconds, as 
is evident in Fig. [3fc. 

4.7.2. Region II - Shocked Wind 

Region II is the shocked wind that has piled up behind 
the shell, where the density jumps by about an order of 
magnitude and transitions from an r~ 2 density profile 
to a nearly flat one (Fig. [3£i). Although the increased 
density has little effect on the injection break frequency, 
both the cooling break and the absorption break each 
change abruptly by approximately an order of magnitude 
(Fig. Et)- There is a moderate change in the flux at all 
frequencies, though the magnitude of the change at each 
frequency and whether it is an increase or a decrease de- 
pends upon the values of the break frequencies. Because 
region IPs imprint on the light curve is modest at best 
even when the emitting region's spherical nature is ig- 
nored, it is likely that the inclusion of spherical emission 
would eliminate the imprint entirely. This is because the 
jet is still highly relativistic when it enters region II, and 
the observer does not see the radiation that is emitted 
when the jet enters region II all at once. The result is 
that the small bump in the light curve that occurs at the 
interface between regions I and II would be smoothed 
out, and region II would have essentially no effect on the 
light curve. 

Depending on the progenitor's mass-loss rate, the mass 
of the shell, and the delay between the expulsion of the 
shell and the burst, the jet can arrive at this region in 
several minutes in the limit of low mass-loss rates, high 
shell masses and short delays between the shell ejection 
and the burst to a year in the limit of high loss rates, 
low shell mass and long delays between the shell ejection 
and the burst. In Fig. [4] we show how the time between 
shell ejection and the burst governs when the jet reaches 
region II, and hence when the light curves and breaks 
would diverge from those expected for r~ 2 wind profiles 
or uniform density fields. In Fig. [5^i we show how mass- 
loss rates and shell masses impact these arrival times. 

4.7.3. Region III - Dense Shell 

Region III is the dense shell ejected by the progeni- 
tor. The jet collides with it in about an hour to several 
years after the burst, depending upon the progenitor's 
wind mass-loss rate, the mass of the shell, and the de- 
lay between shell ejection and the burst (Figs. |4] & [5^)- 
When the jet crosses into region III it abruptly becomes 
non-relativistic because of the large density jump there, 
which can be up to ten orders in magnitude (Fig. [SJl) . A 
bright, highly relativistic reverse shock almost certainly 
forms, which we have neglected for simplicity. 

The transit time through the shell depends on its mass 
and the time between its ejection and the burst. Mas- 
sive shells decelerate the jet more than less massive ones, 
increasing the time that the jet is within the shell. The 
time between the ejection and burst governs the degree 
to which radiative cooling flattens the shell into a thin, 
cold, dense structure, shortening the time the jet is inside 
the shell. Both factors cause the crossing time to vary 
from as little as a day to as much as several hundred days 
(Fig. [Sh,). 

The shell leaves a clear imprint on the spectrum of the 
jet. Upon collision, there is a sharp drop in the cool- 
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ing and injection break frequencies and an increase in 
the absorption break frequency. This causes an abrupt 
increase in the flux at all frequencies (Fig. [SJd). Inside 
the shell, the cooling break frequency rises and the mag- 
netic field strength falls as the jet decelerates, and the 
injection break evolves to lower frequencies. 

The delay in arrival time of photons emitted away from 
the observer's line of sight will certainly decrease the im- 
print of the circumstellar shell on the light curve. At x- 
ray frequencies and higher, where the imprint of the shell 
is least dramatic, there may be little if any increase in the 
flux at the time that the jet enters the shell, though there 
should be a break in the light curve at the time that the 
jet enters the shell where the flux will begin to decay at a 
slower rate. The jet becomes nonrelativistic soon after it 
enters the circumstellar shell, greatly diminishing the im- 
portance of the spherical nature of the emitting region. 
As the jet traverses the shell, the observed flux must 
therefore approach the value that it would have if spher- 
ical emission was ignored. The maximum flux F v occurs 
at some frequency t^fi = min(f m , v c ) in the circumstel- 
lar shell if v a < min(V m , z/ c ), and at ^"hcfi = v a otherwise. 
The increase in flux is largest for frequencies below i^K, 
(Fig. [S|), and the shell's imprint on the lightcurvc is also 
consequently the largest. Even including the effect of 
spherical emission, a significant rcbrightcning should still 
occur for frequencies below i^jtji (typically optical or in- 
frared frequencies and below). While spherical emission 
will have a strong effect on the lightcurve at the time 
that the jet enters the shell, it will not significantly al- 
ter the light curve when the jet exits the shell because 
the jet is non-relativistic. The abrupt drop in the flux 
that occurs when the jet exits the shell should therefore 
remain clearly visible in the light curve at all frequencies. 

4.7.4. Region IV - Low-Density Cavity 

Region IV is the extremely low-density cavity created 
when the fast wind beyond the shell detaches from and 
races ahead of it. The densities in this region are so low 



(io- 



10 8 cm 3 , Fig. [Sh) that the flux density drops 



at all wavelengths when the jet exits the shell (Figs. 
&[5}d). By this stage, the jet has become non-relativistic, 
and photons emitted away from the line of sight are no 
longer significantly delayed with respect to those emitted 
by the portion of the jet that moves directly toward the 
observer. The observer sees the entire leading edge of the 
jet reach region IV at nearly the same time, resulting in 
an abrupt drop in the flux of several orders of magnitude 
at all but radio wavelengths, where the drop is not quite 
as large. Because the jet does not sweep up much mate- 
rial in this region, its magnetic field strength and velocity 
taper off slowly. Consequently, the break frequencies are 
roughly constant (Fig. |3Jd) with v m <C v a <C v c , and 
the spectrum is fairly constant for tens to hundreds of 
days (Figs. [3b & [5b). The absorption break generally 
lies below 1 GHz so the flux density becomes negligible 
for frequencies above the radio band (Fig. [3}d). 

4.7.5. Region V - Detached Wind 

Eventually, the jet crosses the rarefied region and 
catches up to the wind that preceded the ejection of the 
shell. The detached wind, region V, exhibits a sudden 
density jump of several orders of magnitude followed by 



an r~ 2 dropoff thereafter (black plot of Fig.[3Ji). Conse- 
quently, another reverse shock may form at the interface 
between regions IV and V. Unlike the reverse shocks ex- 
pected between regions I and II and between regions II 
and III, this one is Newtonian. We evolved the wind and 
the shell out to a radius of 0.3 pc, and only when the 
burst occurs within 100 years of shell ejection does the 
jet overtake the detached wind before it exits the grid. 
Our models predict that the jet will not reach region V 
for at least several years after the burst, which is why we 
do not show it in Fig. [3ji. 




10 10 10 
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Fig. 4. — Density encountered by the jet as it passes through the 
shell as a function of observer time t f, s for shells ejected 100, 200, 
and 500 yr prior to the GRB (left to right). Here, the progenitor 
had a mass loss rate of 10 -6 Mq yr -1 and a hydrogen shell of 
mass 0.1 Mq. These plots, and our calculations, implicitly include 
the slowing of the jet in the dense shell. 

5. DISCUSSION AND CONCLUSIONS 

Our calculations show that GRB light curves in dense 
shells exhibit clear departures from those in canonical 
winds and uniform densities. These features can broadly 
discriminate between classes of GRB collapse scenarios. 
Analytical models predict that light curves for bursts in 
uniform or r~ 2 density profiles are piecewise power law 
segments separated by breaks. In contrast, light curves 
for GRBs in dense shells ejected by the progenitor (pre- 
dicted for both single and binary He mergers) initially 
follow those for simpler environments, but deviate from 
them in most cases on timescales of a few hours to a few 
days. The first departure is a sudden change in the flux 
of about an order of magnitude when the jet enters the 
shocked wind piled up behind the dense shell. A second, 
much more significant departure occurs soon thereafter 
when the jet collides with the dense shell. 

The features discussed above are distinct signatures of 
a GRB in the dense shell associated with some collapse 
scenarios. The GRB 091127 lightcurve may provide ev- 
idence for a jet en countering a thick shell of material. 
iFilgas et al.l ( 201 ID determined that the temporal evolu- 
tion of the GRB 091127 cooling break frequency could 
only be consistent with the standard fireball model if the 
GRB occured in an r 11 medium. Such a sharp increase in 
density with radius is easily produced in a wind bubble 
environment at the trailing edge of the shell. 
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Our work also provides an alternative mechanism 
for the bright flares that arc sometimes seen in GRB 
lightcurves within the first few thousand seconds of the 
onset of the afterglow emission. GRB 081029 produced 
a flare in the optical band that began at ~ 0.035 days, 
peaked at about 0.07 days after a total increase in bright- 
ness of 1.1 magnitudes, and th en slowly dimmed until 
~ 0.2 days (|Nardini et al.ll2011[ ). The structure and du- 
ration of the GRB 081029 optical flare is somewhat sim- 
ilar to the feature in Fig. [3t> that peaks at 0.3 days. 
A less-intense flare was observed in the GRB 071112C 
X-ray afterglow. The fact that no corres ponding flare 
was d etected in the optical was used by iHuang et al.l 
(|2012f) to argue that different mechanisms were respon- 
sible for the time evolution of the X-ray and optical af- 
terglows. The simple power law structure of the optical 
light curve was shown to be consistent with the external 
shock model, whereas late injection and internal shocks 
inside the jet itself were invoked to explain the excess X- 
ray emission. Our work shows that there can be promi- 
nent, chromatic features in the light curves due solely 
to variations in the circumburst medium. For example, 
Fig. 3b clearly illustrates that the structure of the light 
curves vary with energy band, and that it is not neces- 
sary to invoke delayed energy injection or other effects to 
account for these differences. A similar case is that of the 
"giant" X-ray flare of GRB 050502B, wherein the pho- 
ton count rate was seen to in crease by a factor of 500 at 
345±30s (jFalcone et al.ll200l . Like GRB 071112C, GRB 
050502B was not ob served to have an optical counterpart 
(|Falcone et al.ll2006h . which has been used as evidence of 
a late-injection event. Our models can produce flares of 
comparable magnitude and duration, but can also ex- 
plain the broad plateau in the X-ray flux that lasts from 
the end of the flare up until about a day after the burst. 
A GRB jet that encounters a thick shell of circumburst 
material will produce a flare that will rapidly dim as the 
jet becomes nonrelativistic. The flux then remains nearly 
constant until the jet emerges from the leading edge of 
the roughly uniform-density shell and encounters a sharp 
drop in the medium density. 

The unique imprint of dense shells on GRB light curves 
can be used to constrain their properties. As shown in 
Fig. [SJd, the afterglow flux is somewhat sensitive to the 
ambient density out to ~ a day for a shell ejected 500 
yr before the GRB. The light curve is also sensitive to 
the mass of the shell. The crossing time in a 0.1 M Q 
shell is a day or less, after which there is a sharp drop in 
the flux. On the other hand, a 1.0 M Q shell has crossing 
times of tens to hundreds of days, which creates a broad 
plateau in the light curve. The mass-loss rate of the 
progenitor and the time at which the jet reaches the shell 
are also manifest in the gamma-ray light curve (Fig.[5jD). 
These collectively constrain the mass of the shell and the 
properties of the wind before and after the ejection. 

It is difficult to distinguish exact GRB progenitors be- 
cause many progenitors predict very similar mass ejec- 
tions. The morphology and location of the dense shell 
at the time of the burst are determined by three fac- 
tors: the delay time between shell ejection and the burst, 
the mass of the shell, and the wind mass-loss rate of 
the progenitor. In general, binary mass ejecta will be 
slower, but more massive, than stellar eruptions. Tim- 
ing alone of a flare, however, does not provide a unique 
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Fig. 5. — Density profiles (a) and gamma-ray light curves (b) 
(3x 10 8 GHz, or 1.24 keV) for a GRB occurring 500 yr after ejection 
of the progenitor's hydrogen shell. Black: rn w = 10 -6 Mq yr — 1 ; 
blue: rh w = 10 -5 M yr -1 ; red: rh w = 10 -4 M© yr -1 . Solid 
lines: 0.1 Mq shells; dotted lines: 1 Mq shells. 

constraint (the position of the shell is a function of its 
velocity and time between ejection and collapse). How- 
ever, some progenitors predict specific structures. The 
helium-merger model, for example, predicts a massive 
shell very close to the exploding star. The first possible 
evidence of s uch a progenitor m ay be the recent "Christ- 
mas" burst (jThone et al.ll2011h . With detailed models, 
we may be able to place velocity and mass constraints 
on these shells. With such information, we will both 
be able to better understand massive star evolution and 
constrain the progenitors of GRBs. 

In our models we have adopted some approximations 
and neglected some effects on our light curves. First, 
we neglect the spherical nature of the emitting region, 
which, if included, would un doubtedly result in a modi- 
fication of our light curves (|Fenimore et alj [1996) . We 
have investigated this effect, and the dominant result 
is that the light curve is broadened and becomes more 
smooth, decreasing the imprint of the circumstcllar shell. 
For the majority of the duration of the jet's passage 
through the shell, however, the jet is non-relativistic, 
which will tend to diminish the effect. The shell's im- 
print on the light curve, though lessened, remains de- 
tectable, especially at frequencies where the imprint of 
the shell is largest when spherical emission is ignored. 
Reverse shocks might form at the interfaces between re- 
gions I and II, regions II and III, and regions IV and V, 
and could seriously aff ect the flux of the burst (see i.e. 
iNakar fc Granotl 120071 ) at those radii. Since magnetic 
fields in GRB jets are not well understood, it is also a 
simplification to assume that they are in equipartition. 
Finally, we present synchrotron light curves only. Inverse 
Compton scattering may become important after the jet 
emerges from the dense shell and could increase the flux 
at high frequencies and late times. 

That said, our semi-analytical method can broadly dis- 
criminate between progenitors of GRBs and compute ap- 
proximate GRB light curves and light curves in general 
density fields. Consequently, it is applicable to many 
other GRB host environments besides uniform media 
and winds. Our modeling of the circumburst environ- 
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ment with the Zeus-MP code shows a wide departure 
from the canonical power law density models, with sud- 
den jumps in density of up to ten orders of magnitude 
that have a measureable effect on the light curve. Efforts 
are now underway to apply this method to model obser- 
vational signatures and detection thresholds for Popu- 
lation III ga mma-ray bursts in primordial H II regions 
at z ~ 20 dWhalen et all 120041: iKitavama et all 120041; 
lAlvarez et al.ll2006HAbel et alj|2007t iWise fc Abe1l2008f) 
and both P op II and Pop III GRBs in primeval galaxies 
at z ~ 10 (jWise et al.l 120121 ) for potential successors to 
Swift, suc h as the Join t Astrophysics Nascen t Satellite 
(JANUS, IRomind [20081 IBurrows et all f2010h and Lob- 
ster. 

Our method provides a means to obtain a qualitative 
understanding of gamma-ray burst light curves. For a 



more precise calculation, we must turn to computer sim- 
ulations. Both special-relativistic magnetohydrodynam- 
ical simulations and particle-in-ccll (PIC) calculations of 
GRB jets in circumburst media are now under devel- 
opment, and will soon reveal the light curves of these 
cosmological explosions in unprecedented detail. 
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